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An updated determination of the parameters of the Cabibbo-Kobayashi-Maskawa matrix is presented. 


1. Introduction 

In the Standard Model, weak interactions 
of quarks are governed by the four param¬ 
eters of the CKM matrix [|^ which, in the 
Wolfenstein-Buras parametrisation are la¬ 

belled as A, A, p and fj. Measurements of 
semileptonic decays of strange and beauty par¬ 
ticles are the main sources of information on A 
and A, respectively. The values of |, |I4f,/I4fe|, 
Arrid and Arris provide a set of four constraints 
for p and fj. These constraints depend, in addi¬ 
tion, on other quantities obtained from measure¬ 
ments and/or theoretical calculations. The re¬ 
gions of p and fj preferred by the four constraints 
are expected to overlap, as long as the Standard 
Model gives an overall description of the vari¬ 
ous experimental observations. In this paper, we 
summarize the results of our analysis of the CKM 
matrix. Further details on various aspects of this 
analysis can be found in ref. [Q. 


can be related to the mass difference be¬ 
tween the light and heavy mass eigenstates 
of the system 

Anid = VcSixt) A^\^ [( 1 - 

P?+fj^] , ( 2 ) 

where S{xt) is the Inami-Lim function j^] 
and Xt = rri^/M^. rrit is the MS top mass, 
and pc is the perturbative 
QCD short-distance NLO correction. The 
remaining factor, Jb^^Bb^, encodes the in¬ 
formation of non-perturbative QCD. Apart 
for p and fj, the most uncertain parameter 

in this expression is \JBb^- The value 
of Pc = 0.55 ± 0.01 has been obtained in ||] 
and we used mt = (167 ± 5) GeV, as de¬ 
duced from measurements of the mass by 
CDF and DO Collaborations 0. 


2. Basic Formulae 

Four measurements restrict, at present, the 
possible range of variations of the p and fj pa¬ 
rameters: 


• The limit on the lower value for the time 
oscillation period of the B^ — 5° system is 
transformed into a limit on Anig and com¬ 
pared with Arud 


• The relative rate of charmed and charmless 
5-hadron semileptonic decays which allows 
to measure the ratio 


Vcb 


1 - ^ 
^ 2 


y/p 


/i2 _i_ ^2 


( 1 ) 


• The B^ — B^ time oscillation period which 


Arud _ ruBaf^^BBa f A \ ^ 
Am, “ tubJIBb^ Vl-^/ 

[{l-pf + f]. (3) 


The ratio ^ = Jb, \JBbJ /b, \JBb^ is ex¬ 
pected to be better determined from the¬ 
ory than the individual quantities entering 
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into its expression. In our analysis, we ac¬ 
counted for the correlation due to the ap¬ 
pearance of Arrid in both Equations (|^) and 

(i- 


• CP violation in the kaon system which is 
expressed by \£k\ 


kifi 


Ce fj -riiS{xc) + 

mSixt) {A^\^ (1 - p)) + 

r]3S{xc,Xt) Br , 


(4) 


where 


Atur 


S{xi) and S{xi,Xj) are the appropriate 
Inami-Lim functions § of Xq = 
including the next-to-leading order QCD 
corrections The most uncertain pa¬ 

rameter is Br. 

Constraints are obtained by comparing present 
measurements with theoretical expectations us¬ 
ing the expressions given above and taking into 
account the different sources of uncertainties. In 
addition to p and fj, these expressions depend on 
other quantities which have been listed in Table 
1^. Additional measurements or theoretical deter¬ 
minations have been used to provide information 
on the values of these parameters. 


3. Inferential framework 

The phenomenological analysis is performed 
using the Bayesian inference. In this framework, 
every parameter entering the constraints, regard¬ 
less their theoretical or experimental origin, is 
characterized by a probability density function 
(p.d.f.). We assign the p.d.f.s of the different pa¬ 
rameters as shown in Table |^. These distributions 
are taken as Gaussian or flat or a convolution of 
the two, according to the origin of the uncertainty 
being purely statistical or coming from influence 
quantities (such as theoretical parameters or sys¬ 
tematic errors in experiments). The parameters 


p and fj have also an a-priori p.d.f. which is as¬ 
sumed to be flat. Using these distributions and 
the experimental constraints discussed in the pre¬ 
vious section, we can build an overall likelihood 
and obtain a-posteriori p.d.f.s for p and 77 or any 
other quantities of interest. 

This method provides a theoretically-sound 
approach that allows a consistent treatment of 
the systematic and theoretical uncertainties and 
makes it possible to define regions where the val¬ 
ues of p and fj (as well as of other quantities) are 
contained with any given level of probability. 

The Bayesian method applied to the CKM- 
matrix analysis is discussed at length in ref. 0 . 

A different analysis based on frequentistic tech¬ 
niques can be found in ref. [T^ . This paper 
also puts forth a rather academic argument that, 
in the mind of its authors, should demonstrate 
that the Bayesian method unfairly narrows the 
region of predicted results for quantities depend¬ 
ing on more than one theoretical parameter. The 
argument is the following: consider a quantity 
Tp = X 1 X 2 X 3 ... xn, defined as the product of N 
“theoretical” parameters Xi, and assume to know 
that these parameters lie in the range [—1,1]. Fol¬ 
lowing the Bayesian method one easily finds that 
the resulting p.d.f. of Tp peaks at zero and be¬ 
comes more and more peaked as the number N of 
parameters increases. This leads the authors of 
refs. I^Jl^ to conclude that the Bayesian method 
is “dangerous”, since a safe method should have 
predicted Tp in the range [—1,1]. However, our 
conclusion is quite the opposite. It is perfectly 
natural that the p.d.f. of Tp peaks as N increases, 
being simply the effect of combinatorics (unless, 
of course, there are reasons to believe that the 
different determinations of the Xi are correlated). 
The singularity of the p.d.f. as N goes to infin¬ 
ity, which has been pointed out as a pathology 
of the Bayesian method, has also a simple expla¬ 
nation: the knowledge of an infinite number of 
p.d.f. for the Xi corresponds to determining Tp 
with infinite precision. Therefore we believe that 
the argument presented in ref. II does not affect 
the validity of the Bayesian method in any way. 
We rather find that methods which predict Tp in 
the range [— 1 , 1 ] are deliberately throwing away 
information. This may be justified and reason- 
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Table 1 

Values of the quantities entering into the expressions of \£k\i and Arris. In the third and 

fourth columns the Gaussian and the flat part of the uncertainty are given, respectively. 


Parameter 

Value 

Gaussian 

a 

Uniform 

half-width 

Ref. 

A 

0.2237 

0.0033 


4 

1 

|Vcb| 

40.7 X 10-3 

1.9 X 10-3 


4 

1 

\Vub\ 

36.1 X 10“^ 

4.6 X 10-^ 

- 


4 

1 

\£k\ 

2.271 X 10-3 

0.017 X 10-3 

- 

J 

E 

1 

Arrid 

0.489 ps-3 

0.008 ps-3 

- 

1 

1C 

1 

Arris 

> 14.6 ps-3 at 95% C.L. 

see text 


1C 

1 

mt 

167 GeV 

5 GeV 

- 


7 


mb 

4.23 GeV 

0.07 GeV 

- 

m 


me 

1.3 GeV 

0.1 GeV 

- 


E 

1 

Bk 

0.87 

0.06 

0.13 


4 


fsd\j Bsi 

230 MeV 

25 MeV 

20 MeV 


I 

C _ SBs s/Bbi, 

1.14 

0.04 

0.05 


I 

fSd 






as 

0.119 

0.003 

- 


s 

1 

m 

1.38 

0.53 

- 


s 

1 

V2 

0.574 

0.004 

- 


C 

1 

V3 

0.47 

0.04 

- 


s 

1 

Vb 

0.55 

0.01 

- 


c 

1 

Ik 

0.159 GeV 

fixed 


E 

1 

Amx 

0.5301 X 10-2 ps-3 

fixed 


E 

1 

Gf 

1.16639 X 10-3 GeV-2 

fixed 


E 

1 

mw 

80.42 GeV 

fixed 


E 

1 

mBO 

5.2792 GeV 

fixed 


E 

1 

rriBO 

5.3693 GeV 

fixed 


E 

1 

mx 

0.493677 GeV 

fixed 


9 
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able in specific and well-motivated cases, but as a 
general rule, we find it quite contrary to the spirit 
of this kind of analyses aiming to use present the¬ 
oretical and experimental information to extract 
the best determination of the unitarity triangle in 
the Standard Model. 

4. Input Parameters 

The values of the input parameters used in the 
anaysis are collected in Table |^. For all the the¬ 
oretical parameters, we have used results taken 
from lattice QCD. There are several reasons for 
this choice, which has been adopted also in previ¬ 
ous studies of the unitarity triangle IMl- Lat¬ 


tice QCD is not a model, as the quark model for 
example, and therefore physical quantities can be 
computed from first principles without arbitrary 
assumptions. It provides a method for predicting 
all physical quantities (decay constants, weak am¬ 
plitudes, form factors) within a unique, coherent 
theoretical framework. For many quantities the 
statistical errors have been reduced to the percent 
level. Although most of the results are affected 
by systematic effects, the latter can be “system¬ 
atically” studied and eventually corrected. All 
the recent literature on lattice calculations is in¬ 
deed focused on discussions of the systematic er¬ 
rors and studies intended to reduce these sources 
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Figure 1. The likelihood ratio R(Arris) used in 
the analysis. 


of uncertainty. Finally, in cases where predictions 
from lattice QCD have been compared with ex¬ 
periments, for example fus i the agreement has 
been found very good. Obviously, for some quan¬ 
tities the uncertainty from lattice simulations is 
far from being satisfactory and further effort is 
needed to improve the situation. Nevertheless, 
for the reasons mentioned before, we think that 
lattice results and uncertainties are the most re¬ 
liable ones and we have used them in our study. 

Results from the LEP working groups have 
been used for Am-d and Arris. LEP and CLEO 
measurements of \Vcb\ and |t4f,| have been com¬ 
bined to obtain the values in Table ^ Details on 
the methods used for combining the various mea¬ 
surements can be found in Q] . In order to include 
the information from the lower limit on Arris iu 
the analysis, we have found it better to use the 
likelihood ratio R defined as 


R(Ams) = , 

£(oo) 


( 6 ) 


with 




(7) 


where A is the measured oscillation amplitude at 
a given value of Arris, expected to be equal to one 
at the physical Arris and to vanish elsewhere. The 
likelihood ratio is more effective than the usual 
likelihood as it exploits the fact that the oscilla¬ 
tion is signaled by the amplitude being both com¬ 
patible with one and incompatible with zero. A 
thorough discussion of this method is presented 
in Q| . The likelihood ratio R used in the analysis 
is shown in Eigure ^ 


5. Results 

The region in the (p, fj) plane selected by the 
measurements of \eK\, |Luf)/Rcfe|, ^rnd and from 
the information on Arris (using the R function in 
Figure |^) is given in the upper part of Figure |^. 
On the lower-left part, the uncertainty bands for 
the quantities, obtained using Equations ( 0 )-®, 
are presented. Each band, corresponding to only 
one of the constraints, contains 68 % and 95% of 
the events obtained by varying the input parame¬ 
ters. This comparison illustrates the consistency 
of the different constraints provided by the Stan¬ 
dard Model. The measured values of p and fj are 

p = 0.218 ± 0.038, fj = 0.316 ± 0.040 ( 8 ) 

The two quantities are practically uncorrelated 
(correlation coefficient of -5%), as it can be seen 
from the contour plots in the (p, p) plane. Eitted 
values for the angles of the unitarity triangle have 
been also obtained 

sin(2 P) = 0.696 ± 0.068 , 7 = (55.5 ± 6.2)° , 

sin(2a) =-0.42±0.24. (9) 

In Eigure the p.d.f. for the angles of the 
unitarity triangle are given. Few comments can 
be made: 

• sin(2/3) is determined quite accurately. This 
value has to be compared with the world 
average sin(2/3) = 0.79 ± 0.10 ]I^ ; 

• the angle 7 is known within an accuracy of 
about 10%. It has to be stressed that, with 
present measurements, the probability that 


A log£°“(Ams) 
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P 




Figure 2. Contour plots in the {p,r]) plane. 

7 is greater than 90° is only 0.03%. With¬ 
out including the information from Arris, it 
is found that 7 has 4% probability to be 
larger than 90°. 

As four constraints are used to determine the 
values of two parameters, it is possible to relax, 
in turn, one (or more) of these constraints, still 
obtaining significant confidence intervals. An in¬ 
teresting exercise consists in removing the theo¬ 
retical constraint for Bx in the measurement of 
le^l ([Q-[^). The corresponding selected re¬ 
gion in the (p, fj) plane is shown in lower-right 
plot of Figure where the region selected by the 
measurement of \eK\ alone is also drawn. This 
comparison shows that the Standard Model pic¬ 


ture of CP violation in the K system and of B de¬ 
cays and oscillations are consistent. In the same 
figure, we also compare the allowed regions in the 
(p, ff) plane with those selected by the measure¬ 
ment of sin(2/3) using J/ijjKs events. 

Using constraints from 6 -physics alone the fol¬ 
lowing results are obtained 
f , — n OQ4-1-0.050 

0.167 < p < 0.400 at 95% prob. 

sm{2/3)=0.676t°o°ll (10) 

0.430 < sin(2 j3) < 0.820 at 95% prob. 

Another way for illustrating the agreement be¬ 
tween K and B measurements consists in compar¬ 
ing the values of the Bk parameter obtained in 
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Figure 3. Some a-posteriori p.d.f. 


lattice QCD calculations with the value extracted 
from Equation (^, using the values of p and fj se¬ 
lected by ^-physics alone 

^b-phys ^ o.SStouL 

0.65 < < 1.65 at 95% prob. (11) 

Since is not limited from above, for the 
present study, probabilities are normalised as¬ 
suming Bk <5. 

The importance of B^-B^ mixing can be illus¬ 
trated from the p.d.f. of the angle 7 obtained with 
or without including the Arus constraint (respec¬ 
tively light and dark lines in the lower-left plot of 
Figure ||). High values for 7 are excluded at high 
confidence level by the experimental lower limit 
on Atos. 


It is also possible to extract the probability dis¬ 
tribution for Ams, from which one obtains 

ArUs = (16.1 ± 3.2) ps“^, 

9.4 < ArUs < 23.0 ps“^ at 95% prob. (12) 

If the information from the B^ — B^ analyses is 
included, results become 

Ams = (17.lt(,;g) ps“^, 

15.4 < Ams < 20.3 ps“^ at 95% prob. (13) 

These values are in agreement with the recent es¬ 
timate of Ams = 15.8(2.3)(3.3) ps~^, presented 

in H. 

The value of \jB b^ can be obtained by re¬ 
moving the theoretical constraint coming from 
this parameter in B'^-B^ oscillations. Using 
the two other theoretical inputs, Bk and 
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fsd y ^Bd is measured with an accuracy which 
is better than the current evaluation from lattice 
QCD, given in Table |^. We obtain 

/B.V^=(228±12)MeV. (14) 

The present analysis shows that these results are 
in practice very weakly dependent on the exact 

value taken for the uncertainty on fsd \JBsd- An 
evaluation of this effect has been already pre¬ 
sented in 1 ^^ where the flat part of the theoretical 

uncertainties on fsd\jBsd '^^-s multiplied by two. 
Similar tests will be shown in Section |^. 

It is even possible to remove the theoretical 

constraints on both fsd \JBsd and Bk and ob¬ 
tain the simultaneous lower bound 

Bk>Q.^ and 150 MeV (15) 

at 95% probability. 

6 . Stability of the results 

The sensitivity of present results on the as¬ 
sumed probability distributions attached to the 
input parameters was studied. The comparison 
of the results obtained by varying the size of the 
theoretical uncertainties has been done to eval¬ 
uate the sensitivity to these variations of uncer¬ 
tainties quoted on fitted values. This must not 
be taken as a proposal to inflate the uncertainties 
obtained in the present analysis. 

In these tests, all values for uncertainties of the¬ 
oretical origin have been, in turn, multiplied by 
two. For the quantities \Vub\ and |14b|, new p.d.f. 
have been determined, following the prescriptions 
mentioned in Section ^ and used in the analy¬ 
sis. The main conclusion of this exercise is that, 
even in the case where all theoretical uncertainties 
are doubled, the unitarity triangle parameters are 
determined with an uncertainty which increases 
only by about 1.5. 

7. CKM angle 7 from B Ktt, tttt decays 

It has been recently suggested that it is pos¬ 
sible either to extract the CKM angle 7 using 
the measured BRs of B — s- Ktt, tttt decays or, at 



Figure 4. Dependence of 7 on the maximum 
value of the phenomenological parameter \pa\ ap¬ 
pearing in the model used to compute BR{B —> 
Ktt, tttt). The bands correspond to Icr ranges. 


least, use this experimental information to im¬ 
prove the results of the unitarity triangle anal¬ 
ysis 1^. Although in principle this is certainly 
true, in practice it requires a very good control 
on the hadronic uncertainties entering the the¬ 
oretical predictions of these BRs. Contrary to 
recent claims, however, the theoretical progresses 
in understanding in the infinite mass limit factor¬ 
ization of hadronic matrix elements do not help 
in the specific case. In fact, B Ktt,tttt chan¬ 
nels get large, if not dominant, contributions from 
power-suppressed terms, for which a theory has 
not been developed so far. For this reason, the ap¬ 
proach of ref. 1 ^ , where power-suppressed con¬ 
tributions are factorized together with the leading 
terms, should be regarded as a phenomenological 
model, and one with quite specific assumptions. 
As a tool to extract fundamental parameters such 
as the CKM angles, it has the same difficulties 
in assessing the theoretical uncertainties as any 
other model and therefore should be used with 
caution. 

To further illustrate this point, we have in¬ 
cluded in our analysis the constraint coming from 





















BR{B Ktt, tttt) using the model of ref. ||2^ 
and considered the dependence of the predicted 
value of 7 on the phenomenological parameter pA ■ 
This free parameter is introduced in the model 
to account to some extent for the infrared di¬ 
vergences appearing in the perturbative calcula¬ 
tion of power-suppressed corrections. According 
to ref. 1 ^, the allowed range of variation of pA 
is \pa\ < 1. We find instead that the data pre¬ 
fer larger values \pa\ ^ 4-6. Figure || shows that 
the value of 7 , which we obtain by including con¬ 
straints from non-leptonic i?i?s, deviates from the 
one of eq. (||) only for value of max(|p^|) ~ 1-3. 
On the other hand, when pA is allowed to get 
the values preferred by the data, no new infor¬ 
mations on 7 are obtained besides those coming 
from more reliable constraints. In other words, 
it is just the restricted range for pA adopted in 
ref. [ p5[ which originates the claimed “improve¬ 
ment” in the determination of the CKM parame¬ 
ters due to non-leptonic BRs. To our knowledge, 
this choice has no compelling theoretical or phe¬ 
nomenological basis. 

We rather believe that, at present, our control 
of the hadronic uncertainties in the B Ktt, tttt 
BRs is not developed enough to allow using these 
modes for costraining the CKM parameters. 
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